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Abstract 

In this article, we report on the curious phenomena of anomalous spreading in a system of coupled 
Fisher-KPP equations. When a single parameter is set to zero, the system consists of two uncoupled 
Fisher-KPP equations which give rise to traveling fronts propagating with the unique, minimal KPP 
speed. When the coupling parameter is nonzero various behaviors can be observed. Anomalous spreading 
occurs when one component of the system spreads at a speed significantly faster in the coupled system 
than it does in isolation, while the speed of the second component remains unchanged. We study these 
anomalous spreading speeds and show that they arise due to poles of the pointwise Green's function 
corresponding to the linearizion about the unstable homogeneous state. These poles lead to anomalous 
spreading in the linearized system and come in two varieties - one that persists and leads to anomalous 
spreading for the nonlinear system and one that does not. We describe the mechanisms leading to these 
two behaviors and prove that one class of poles are irrelevant as far as nonlinear wavespeed selection is 
concerned. Finally, we show that the same mechanism can give rise to anomalous spreading even when 
the slower component does not spread. 

MSC numbers: 35C07, 35K57, 34A26 
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1 Introduction 

In this article, we study anomalous spreading speeds in a system of coupled Fisher-KPP equations, 

Ut = du xx + a (u — u 2 ) + fiv(l — u) 

v t = v xx + (v-v 2 ). (1.1) 

We are interested in the following problem. Consider a positive, compactly supported perturbation of the 
unstable homogeneous state (u,v) = (0,0). What is the asymptotic speed of propagation associated to the 
u component? 

Consider first the case when j3 — or v is identically zero. Then the equation governing the dynamics of 
the u component is the scalar Fisher-KPP equation, 

u t = du xx + a(u - u 2 ). 

This equation has been studied in great detail by a number of authors, see [H [TBI UJ H] among others. It is 
well known that compactly supported, positive perturbations of the zero state evolve into a pair of counter 
propagating fronts. These fronts travel with asymptotic speed 2\fda. At the same time, observe that the 
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v component decouples and its evolution is also given by a Fisher-KPP equation with a selected spreading 
speed of two. For /3 > the evolution of the u component depends on that of the v component. One might 
conjecture that the selected speed of the u component for the full system will be either 2 or 2yda, whichever 
is larger. This turns out to not be the case, and there are large swaths of the (d, a) parameter plane for which 
speeds of propagation are observed that exceed both 2 and 2y/da. See Figure |l] for an illustration of those 
parameter values leading to faster speeds of propagation in the u component. This phenomena was first 
observed in [33] and given the label of anomalous spreading. In this paper, we give a complete description 



of the anomalous spreading that arises in system (1.1) and elucidate the mechanisms leading to these faster 
than expected speeds of propagation. 



A natural starting point for the analysis is the linearization of (1.1) about the homogeneous unstable 
zero state. The fastest speed of propagation associated to the linearization is called the linear spreading 
speed and is of interest for many reasons. Practically speaking, the linear spreading speed can be inferred 
from the singularities of the pointwise Green's function. We delay a precise description of this process until 
section [2] but we remark that this was first observed within the plasma physics community; see [5] [3] for 
descriptions of the original work aimed at differentiating between absolute and convective instabilities. From 
a physical standpoint, one would like to understand the relationship between the linear spreading speed and 
the nonlinear spreading speed. This is a complicated issue and a complete overarching theory is not known. 
It has been shown in [TH] HH] that a large class of invasion problems are linearly determinate; that is the 
nonlinear selected speed is exactly the linear spreading speed. These linear determinate fronts are pulled 
fronts according to the commonly used terminology reviewed for example in |23j . This, of course, does not 
hold in general and there exists examples for which the selected wavespeed is not the linear one, see for 
example Nagumo's equation [TT]. These nonlinearly determined fronts are referred to as pushed fronts. 

Another natural question is whether the linear spreading speed places a lower bound on the speed 
of propagation for the nonlinear system. Wavespeed selection for general systems of equations is usually 
described in terms of marginal stability of the invasion front, j6]. The selected front should be marginally 
stable against compactly supported perturbations, which is to say that perturbations of the traveling front 
should neither grow nor decay when viewed pointwise in a traveling frame moving with the speed of that front. 
Since the linear spreading speed describes the spreading speed of perturbations of the unstable homogeneous 
state, one might imagine that a perturbation placed sufficiently far ahead of the front interface should spread 
with the linear spreading speed. This seems to suggest that the linear spreading speed should place a lower 
bound on the speed of propagation for the nonlinear system. This is the case for scalar problems, but 
for systems this bound no longer holds. This fact was recently shown in an example of a Lotka-Volterra 
competition model, [T5] . 

Ut = c u xx + (1 — u — a\v)u 

v t = v xx + r(l - a 2 u - v )v. (1.2) 
Here < e <C 1, r > and a\ < 1 < a 2 . The unstable steady state is the point (u,v) — (0,1). As 



is the case in (1.1), the linearization about this fixed point has a skew-product structure. On the level 



of the linearization, the u component decouples and is described by a Fisher-KPP equation that spreads 
with rate 2eyT — ai. The linear coupling to the v component occurs as an inhomogeneous term in the v 
equation and this together with the fast diffusion of the v component leads to 0(1) spreading speeds for 
this system. A crucial feature leading to this behavior is the unbounded exponential growth in time of the u 
component in the linearized system. For the non-linear system, no such exponential growth is possible and 
so the anomalous linear spreading speed does not persist for the nonlinear system. This should be contrasted 



with the dynamics of system (1.1) as well as the example provided in [53]. In both those cases, anomalous 
spreading speeds are observed in both the linear and nonlinear systems. 

Mathematically, we will show that these anomalous linear spreading speeds arise due to poles of the 
pointwise Green's function. This is a crucial difference between the scalar and non-scalar cases. Systems of 
reaction-diffusion equations can support poles of the pointwise Green's function whereas these singularities 
in the scalar case are always accompanied by a loss of analyticity. Poles of the pointwise Green's function 
are important because they represent the linear spreading of one component that is induced by the coupling 
to a component spreading with a different speed. Based upon the results depicted in Figure [l] as well as the 
examples in |24|,I13] it seems that these poles come in two varieties: one that induces faster spreading speeds 
in the nonlinear system and one that does not. In section [3j we will come to this question and describe 
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two different mechanisms leading to anomalous linear spreading speeds. One such mechanism persists in 
the nonlinear system, while the second type does not. We remark that for poles of the pointwise Green's 
function to exist, the linearization must possess a skew-product structure where some of the components 
decouple from the others. Problems for which this skew-product structure exist arise in many areas of pattern 
formation, ecology and chemistry among others, see [T^l H7J [2U] for recent examples. 

One interesting aspect of this invasion process is the range over which the slower v front has an influence 
over the faster u front. For large values of t, the u component will develop into a traveling front propagating 
in a medium for which v is essentially zero. In fact, viewed pointwise in a frame traveling with the asymptotic 
speed of the u component the v component is converging pointwise to zero as t — > oo. Nonetheless, the v 
component, while small, is non-zero and has certain spatial decay rates. These decay rates may be adopted 
by the u component. In contrast to the mechanism that leads to faster linear spreading speeds in (1.2 1, 



this mechanism is not destroyed by the introduction of a nonlinearity and faster speeds of propagation can 
be observed in the nonlinear regime. This picture is somewhat reminiscent of recent work in |12) where a 
similar long range interaction led to faster observed speeds of propagation. An important distinction exists, 
however, since the decay rate of the KPP front is not sufficient to determine the selected speed for the u 
component in (1.1 1. Instead, the complex process by which the compactly supported initial data evolves 



towards the traveling front profile plays a role in the nonlinear speed selection. This process appears to be 
insufficiently understood. 

That anomalous spreading exists despite the fact that the v component is converging pointwise to zero 
in a frame moving with the anomalous spreading speed suggests that this spreading phenomena should not 
be dependent on the instability of the zero state with respect to perturbations of the v component. In fact, 



this is the case. For example, suppose that the equation governing the v component in (1.1 1 was instead 



v t = v xx - v. 

The zero state is stable and any perturbation of that state will relax to zero. However, in the process of this 
relaxation the v component is non-zero and has certain decay rates as x —¥ oo. In an analogous manner to the 
anomalous spreading for ( |1.1[ ), the effect of this spatial decay can give rise to faster speeds of propagation in 
the nonlinear system. We give an explicit example of this in section [4j An important conclusion to be drawn 
here is that one must use caution when reducing a system of reaction-diffusion equations by setting stable 
or slowly propagating components to some constant value. The reduced system have have very different 
dynamics than the original. 
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Figure 1: The selected linear spreading speed for the linear (left) and nonlinear (right) systems. Here / 
denotes a region where the selected speed of the u component is the pulled speed whereas in region II 

the selected speed is 2. The region A is where anomalous speeds are observed. Note the difference between 
the linear and nonlinear cases - the anomalous speeds arising due to fast diffusion of the u component do 
not persist in the nonlinear regime. 



The paper is outlined as follows. In section [2] we compute the linear spreading speed using the pointwise 
Green's function. In section [3j we discuss under what circumstances anomalous spreading in the linearized 
equation leads to anomalous nonlinear spreading. In section [4j we give an example that shows that this 
phenomena arises even when the component that induces the anomalous spreading is pointwise stable. In 
section |5j we contrast numerically observed spreading speeds with the predictions based upon the linearized 
system. We conclude in section [6] with a short discussion. 
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2 Spreading speeds in the linearized system 



In this section, we compute the linear spreading speed associated to the u component in ( 1.1 ). One can think 
of this speed intuitively as the infimum of all speeds for which a moving observer will outrun the instability. 
Put another way, the linear spreading speed is the largest speed for which a transition from absolute to 
convective instability is observed. This transition can be naturally understood in terms of singularities of 



the pointwise Green's function. After linearizing system ( 1.1 ) about the unstable zero state, we will compute 



these singularities explicitly and calculate the linear spreading speed. 

Preliminaries - linear spreading speeds, the pointwise Green's function and pinched double 



roots of the dispersion relation To begin, we linearize (1.1) about the unstable zero state and consider 
its evolution in a moving coordinate frame. Let £ = x — st. Then the linearization in this coordinate frame 
is 

Ut = dv,££ + su^ + au + j3v 

v t — v^ + sv^+v. (2.1) 

Since the zero state is unstable, compactly supported perturbations of this state will grow in norm and 
spread spatially. However, by considering this evolution in a moving frame an important distinction can 
be made as to whether the instability is absolute or convective in this frame. Let us review the notion of 
convective versus absolute instabilities. For absolute instabilities, the instability grows both pointwise as well 
as in norm. When the instability is convective, the instability is transported away from its original location 
and pointwise decay to zero is eventually achieved despite exponential growth in norm. Mathematically, 
the question of how to differentiate between convective and absolute instabilities has a long history in the 
literature, see [5J [31 HH [ZD 123] for the original work within the plasma physics community and subsequent 
extensions. The key observation is that absolute and convective instabilities can be differentiated according 
to whether the pointwise Green's function, G\(£, — y) has singularities in the right half plane. 



We review some basic information about the pointwise Green's function G\(£ — y) for ( 2.1 ). We refer the 



reader to [25 for more information. The pointwise Green's function is computed as follows. After a Laplace 



transform in time, the linear system (2.1 ) reduces to a system of second order ordinary differential equations 
with inhomogeneous right hand side. These equations depend on A in an analytic fashion and for each A the 
corresponding Green's function for this system of equations is denoted Ga(£ — y). That is, we solve 



C u - A \ ( Gl 1 Gf \ _ ( S(y) 



C v - A ; V Gf Gf J v <%) 



(2.2) 



The solution for general initial data can be recovered by convolution. The pointwise Green's function is 
analytic to the right of the essential spectrum since it is the kernel of the analytic resolvent operator. 
Pointwise decay is achieved if G\ can be analytically continued into the left half plane. Of course, analyticity 
properties of Ga(£ — y) depend on the frame of reference, or equivalently on the value of the wavespeed s. 
For s = 0, Ga(£ — y) has a singularity in the right half plane. As s is increased this singularity moves until 
eventually some value of su n is reached so that all singularities of G\ (£ — y) remain in the stable half plane 
for all s > su n . This value of the wavespeed for which the singularity crosses the imaginary axis is the linear 
spreading speed. That is, we have the following definition. 

Definition 1. The linear spreading speed is 

sii n — sup {Ga(£ — y) has a singularity for some A with Re(X) > 0} . 

In turn, singularities of the pointwise Green's function are related to the eigenvalues of the operators 
£u,v ~ ^- Denote these eigenvalues v~ (A, s) and *v(A, s ) where we will henceforth suppress the dependence 
on A and s. These eigenvalues are roots of the characteristic polynomials, dv 2 + sv + a — A and v + sv+ 1 — A. 
The product of these polynomials is the dispersion relation, d{v, A) = d u (p, X)d v (v, A), roots of which relate 
spatial modes e ux to their temporal growth rate e At , 

d(X, v) = {dv 2 + sv + a - \){u 2 + sv + 1 - A). 
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To the right of the essential spectrum, the eigenvalues v+ v lie to the right of the imaginary axis while the 
stable eigenvalues v~ v lie to the left of the imaginary axis. Analytic continuation of G\ depends on the 
ability to analytically continue the stable and unstable subspaces associated to these eigenvalues into the 
essential spectrum. Potential issues arise at values of the spectral parameter A for which an eigenvalue that 
originated on one side of the imaginary axis aligns with an eigenvalue that originated on the opposite side of 
the imaginary axis. These are referred to as pinched double roots of the dispersion relation, that is, d{X, v) 
has a pinched double root at (A*, v*) if 

d(A*,z/)=0, d v d{\*,v*) = 0, Re v ± {\) -> ±00 as Re(X) -> 00. 

We remark that while pinched double roots may prevent analytic continuation of G\, the existence of a 
pinched double root is not a sufficient condition for the existence of a singularity of G\. This is because 
the pointwise Green's function has an analytic continuation if the stable and unstable subspaces can be 
continued analytically and pinched double roots do not imply that these subspaces are not analytic. This is 
easy to observe in our example. If j3 = 0, the equations are decoupled and it is clear that a value of A for 
which — should not lead to a singularity of G\. Nonetheless, for /3 7^ these pinched double roots 
do give rise to singularities of the pointwise Green's function. We will now explicetly compute the pointwise 
Green's function and show that, in some instances, these pinched double roots induce poles of the pointwise 
Green's function. 



Spreading speeds from the dispersion relation Due to the skew-product nature of the linearization 



in (2.1 1, we can compute the roots of the dispersion relation explicitly. We find 



vt(\,s) = ^±^s 2 -4 + 4A. (2.3) 

Pinched double roots of the full dispersion relation occur for all values of A, v and s for which i/+ — v~ v . 
Each pinched double root induces a spreading speed that can be found by selecting s in such a way that the 
value of A at the double root satisfies Re(X) = 0. Thus, we set A = and solve for s. The spreading speeds 
associated to the u and v components in isolation are easily found by setting the terms inside the roots in 



(2.3) to zero. A third spreading speed, the anomalous spreading speed is found by solving for values of s for 
which — with A again set equal to zero. Doing this, we find the following three speeds 

s u ,lin : = 2Vda 



a- 1 l-d 



l-d V a - 1 

Depending on what parameter values, (d, a), are under consideration the selected spreading speed of the u 
component is one of these three speeds. We have the following result. 

Theorem 1. The linear spreading speed for system is 



2 for a < 2 - d 

= J 2\fda for a > 2 d-\ an ^ > I) (2 4) 



y — -7 + \ — i otherwise 

l—a \ a— 1 



These regions are plotted in Figure [T] This result can be established by verifying that the pinched double 
root calculation that leads to the derivation of these speeds actually induces a singularity of the pointwise 
Green's function. To see this, we must calculate G\. 
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Computing the pointwise Green's function We will now compute G\. As is expected, this function 
will have singularities at the values of A for which pinched double roots exist. Furthermore, we will observe 
that these singularities come in two varieties: either G\ is meromorphic in a neighborhood of the singularity 
or it fails to be analytic due to a branch point of the dispersion relation. 



Recall (2.2). The first component, G^ 1 , satisfies the differential equation (C\ — X)G\ 1 — S(y), for which 



it is readily computed that 



f _i e » u (A)«-y) f or £ >y 

G\ 1 (S-y) = < u 7"" + , . (2.5) 

y ' | —L_^«\T (•*)(£-») for £ < y y ' 



In an analagous fashion, we compute 



e f„ Wa-v) for £ > y 



Gf(£- y) = ^ 7"« , ° . (2.6) 

I -^i e ^W(e-») for £ < w y ' 

It is easily verified that G^ 1 has a singularity and loses analyticity at A = a — f-j while G 22 does so at 
A = 1 — Selecting s — s u .un = 2V da places the singularity of Gj^ on the imaginary axis, while selecting 



-j-. kjcict-Liiig, o — o u ,i m — iy fiasco unc oiiiguiaiiuy ui u A 

s = Su,(m = 2 does the same for the v component. Note that the pinched double roots giving rise to the 
anomalous spreading speed have no effect on the analyticity of G^ 1 or G| 2 . In other words, when (3 = 0, 
the spreading speeds are just these linear spreading speeds and the pinched double root giving rise to the 
anomalous speed is not dynamically relevant. 

However, we are interested in the case of (3 ^ 0. In this event, the non-diagonal element G^ 2 is non 
zero and describes the linear behavior of the u component due to coupling with the v component. We will 
see that the pinched double root v~~ = may be relevant here. To compute G^ 2 we consider the linear 
equation 

(£„ - A)Gf + pGf = 0. 

Skipping the specifics of the derivation, the general solution is computed on either half line (£ > y, £ < y) to 
be 

Al4 V) ~\ C _(A)e^>^> - J + — \ *»tMK-v) for £<y ' 
v v u —Vu d u {fv ,A) 

with the initial conditions enforced at £ = y so that G^ 2 is continuously differentiable. This prescribes the 
constants 

c + (A) = 
c_(A) = 

We make several immediate observations. First, we see that G^ 2 loses analyticity whenever any of the 

, 2 2 

eigenvalues v^ v \\) does. To the right of the branch points A = a — |j and A = 1 — ^- we then have that 
G^ 2 is analytic with the possible exception of those A values for which d u {v^ ,A) = 0, or whenever a spatial 
eigenvalue for the v subsystem, i.e. v~ is equal to a spatial eigenvalue for the u subsystem. At any double 
root that does not satisfy the pinching condition, the pointwise Green's function has a removable singularity 
and therefore these points do not pose an obstacle to analytic continuation of G\. On the other hand, 
pinched double roots where = lead to poles of the pointwise Green's function. The location of these 
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poles can be found exactly as was done above, by computing those values of A for which the roots in (2.3) 
form pinched double roots. This establishes Theorem [T] 



3 Anomalous spreading speeds in the nonlinear system 

Whenever (A) = v^(\) we find that the pointwise Green's function has a pole. However, in numerical 
simulations we observe that it is only the poles wherein i/+ = v~ that give rise to faster observed speeds 



G 



of propagation in the nonlinear system, recall figure [T] In this section, we will investigate the physical 
mechanisms that enforce the existence of these poles and explain why these speeds persist to the nonlinear 
system only when the double root involved is z/+ = v~ . In the first case (f+ = v~), we argue that faster 
spreading speeds are observed due to the alteration of the decay rates of the u component through coupling 
to the v component. In the second case (y~ = we argue that these poles arise due to a combination of 
the slower growth but faster diffusion associated to the u component relative to the v component. We will 
also prove that this second type of double root pertains to spreading speeds that are not attainable in the 
nonlinear system. 

The key observation is that the spatial eigenvalues v~ can also be interpreted as giving spatial decay 
rates of some traveling front. These decay rates fall into one of two groups - those with strong decay 
and those with weak decay (v^)- For the nonlinear system, it is known that any front with weak decay is not 
attainable from compactly supported initial data, see for example 23J. As a result, we will show that any 
pinched double roots that involve v£ are not relevant for the nonlinear system. Nonetheless, these double 
roots lead to poles of the pointwise Green's function and hence anomalous linear spreading. By contrast, 
the double roots that involve v~ are still relevant for the nonlinear system since any initial condition with 
this exponential decay rate is dynamically relevant. 

Before discussing these two mechanisms in more detail, we will first recall the interpretation of the 
dispersion relation that relates spatial decay to spreading speeds. Consider the linearized equation, 

Ut = du xx + au. 

As we did when deriving the dispersion relation above, take the ansatz u{t. x) — e xt e ux and substitute this 
into the linear equation. Assuming that v £ M and v < 0, then we find that A = dv 2 + a. From this the 
group velocity, 

A l a /Q 1 \ 

s Q = = dv H — , (3.1) 

v v 

gives the speed with which this exponential propagates. This relationship can be visualized in (y, s) space, 
see Figure 5j For each value of s > 2\fda there are two possible decay rates that propagate with speed s. 
These values of v are denoted f„ and called spatial eigenvalues as they are exactly the eigenvalues of the 
fixed point at the origin for the traveling wave equation, 

du^ + SU£ + u — u 2 = 0. 

Note that the linear spreading speed, 1\fda, is the value of s where these two branches meet. The upper 
branch corresponds to traveling front profiles with weak decay. These fronts are stable in an exponentially 
weighted space, see [22], but are not accessible from compactly supported initial data. On the other hand, 
if the initial data is not compactly supported but instead has an exponential decay rate belonging to this 
weak branch then the initial data will evolve into a traveling front moving with the speed prescribed by the 
group velocity calculation above. For the Fisher-KPP equation, no nonlinear traveling front solutions exist 
that approach the origin with decay given by the branch v~. Therefore, initial data with this strong decay 
evolves towards the KPP traveling front. 

The same general picture holds for the v component. Consequently, we see that pinched double roots 
involving v u and v v coincide with values in the (y, s) plane for which the branches v£ and v~ (or v~ and 
uj) intersect. These anomalous spreading speeds are always observed in the linear system. However, when 
nonlinearities are introduced only the double roots corresponding to v+ = v~ that give rise to anomalous 
spreading speeds in the nonlinear system. We give physical interpretations of the mechanisms leading to 
these anomalous linear speeds and use these interpretations to explain why in one case = v~) they 
persist in the nonlinear system, while in the other (y~ = z/+) they do not. 

• Adoption of weak decay by the u component - In the nonlinear system, the v component is 
evolving towards a traveling front profile moving with speed 2. Associated to this traveling front 
profile is a selected decay rate, here £,e~^. When this decay rate is shallower than the decay naturally 
selected by the u front in isolation (£e~V^) then it is possible that the effect of the linear coupling 
will be sufficient to establish weaker decay rates in the u component than the selected one. These 
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Figure 2: Plots of the group velocity s g as a function of the spatial decay rate v. The solid line is this curve 
for the v component. Note the branch point at the linear spreading speed of 2 with selected spatial decay 
rate of —1. The dashed line is the decay rate - group velocity plot for two different values of (A, nu). On 
the left is the curve for d — 0.1 and a = 5. Note the pinched double root that occurs involving f+ and u~ . 
This pinched double root leads to faster observed speeds of propagation in both the linear and nonlinear 
systems. On the right is the curve for d = 5 and a = 0.1. Here, the system exhibits anomalous spreading 
for the linear system but not for the nonlinear system. 



weaker decay rates lead to faster than expected speeds of propagation. Furthermore, since it is only 
the decay rate that matters, this influence can have an effect over a very large spatial scale. 

We emphasize that the mechanism leading to faster propagation speeds is more subtle than simply 
the adoption of the weak decay rates associated to the v front. Faster speeds of propagation are 
observed only if there exists a double root at z/+ = v~ . The observed spreading speed is then the linear 
spreading speed associated to this double root. In general, this speed is slower than the spreading 
speed associated to weak fronts with decay £e~^. This is due to the rather complicated route that 
compactly supported initial data take in the convergence to their asymptotic form of a traveling front 
solution. 

• Weak growth and fast diffusion in the u equation - There are also double roots and poles of the 
pointwise Green's function associated to double roots wherein v~ = v+ . These poles give rise to faster 
speeds of propagation in the linearized equation, but not in the nonlinear equation - see Theorem [2] 
below for a proof of this in the nonlinear context. 

That anomalous spreading speeds are observed in the linearized equation is apparent from the pointwise 
Green's function. The physical mechanism leading to this phenomena can be best conceptualized in 
the limit of large diffusion in the u component. Here, the v component grows and diffuses as usual. 
These dynamics enter into the u equation as an inhomogeneous source term. The u dynamics adopt 
the exponential growth of the v component, but the diffusion is much larger for the u component and 
this leads to faster linear spreading speeds. Note that this was the same mechanism discovered in [13] 



for the Lotka-Volterra competition model, (1.2) 



Observe that this mechanism depends critically on the exponential in time growth of the v compo- 
nent. For the nonlinear system no such growth is possible and therefore this type of anomalous linear 
spreading does persist for the nonlinear system. This fact is made rigorous in the following Theorem. 



Theorem 2. Consider system (1.1). Suppose that d > 1 and a < 2 d-i - Then the selected wavespeed of the 



u component is bounded above by max{2, 2\/ da} 



Proof: Let s > max{2, 2y/da}. We will construct a front-like super-solution that propagates to the 
right with speed s, see [8]. To begin, it is known that for any C v > 0, 

v(t,x) = mm{l,C v e^ {0 - s){x - st) } 

is a super-solution for the v dynamics. For every x, the u component converges pointwise as t — ¥ oo to the 
KPP front. That being said, the super-solution above places an upper bound on the spatial decay rate at 
any fixed time t > 0. Thus, the evolution of the u component conserves the steepness of the initial data, see 
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[7j. We then turn our attention to the u component. Let 

N[u) = Ut — du xx — au — /3v(l — u) + au 2 . 

We seek u(t, x) so that N(u) > for all i > and i£R. Consider 

-/ \ / 1 for (a; — si) < r 

= I Cue"- ( x " st ) + CiKe"^*-**) for (a; - si) > r ' ( } 

We will show that we may select r and C„ > so that u(t, x) is a super-solution. To begin, one easily verifies 
that u — 1 is a super solution. If t is taken to be larger than r m j n := (i^) -1 logC„ then for x — st > t we 
have that N(u) is given by 



N(u) =u t - du xx ~au- /3C v e-^ {x - st) + 0(v - v) + (3vu + 

Consider 



au 2 . 



du xx + au + (3C v e-^ {x - st) . 



This equation has the solution, 



for 



l(t,x) = C u e"^ x - St) + C v Ke"« ( x ~ st \ 

» = ^ <0. 

d{v v Y + sv v + a 

For (d, a) under consideration, we also have that v~ < v~. As a result, for any C u > we have that 
u(t, x) > for x sufficiently large. Furthermore, u(t, x) has a unique maximum at 

log( 



Kmax Z Z i 

V u - Vv 

for ^ = x — st. For C u > and sufficiently large we have that both £ max > T m i n and the value of u(t, x) 
at this maximum exceeds one. Let r be chosen so that w(i,r) = 1. Then u(t,x) in (3.2) is continuous. It 
remains to verify that N(u(t,x)) > for x — st > r. On this semi-infinite interval we have, 

N(u(t, x)) = N(u(t, x)) = I5{v -v) + /3uv + au 2 > 0, 

establishing u(t, x) as a super-solution. 

This construction works for any s > max{2, 2\J da} and therefore establishes that no anomalous spreading 
speeds are attainable for the nonlinear system in this regime. 



The relevant case An analogous result for the relevant case is more elusive. Suppose that a pinched 
double root exists at — v~ . One might be tempted to take a dynamical systems approach and seek 
traveling front solutions. That is, let Vkpp(x — 2t) denote the Fisher-KPP traveling front solution. The 
dynamics of the v component evolve towards this traveling front solution and it therefore might seem natural 
to study invasion speeds in the inhomogeneous scalar equation, 

u t = du xx + a(u — u 2 ) + I3Vkpp{x — 2i)(l — u). (3-3) 

Taking this approach, however, will over-estimate the actual speed of propagation. The effect of the inho- 



mogeneity in (3.3) is to impart the decay rate of the KPP front onto the u component. The selected speed 



of propagation for (3.3) can then be determined by the group velocity calculation in (3.1). This speed will 
be larger than the speed selected by compactly supported initial data. 

This discrepancy is due to the fact that while for every £ = x — 2t the v dynamics converge to the 
traveling front profiled Vkpp{£,) as t — > oo, the convergence is not uniform. Constructing a sub-solution 
that shows that the nonlinear anomalous spreading speed is always achieved would require some detailed 
knowledge concerning the route that the v component takes during its convergence to the traveling front 
Vkpp{£,)- This does not appear to be completely understood in the literature and is beyond the scope of 
this paper. 
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4 A stable system inducing anomalous spreading 



Viewing the evolution of ( |2.1[ ) in a frame of reference traveling at the anomalous spreading speed, we 
observe that the v component is convectively stable. That is, for any fixed £ = x — s anom t the v component 
converges pointwise exponentially fast in t to the zero state. Taking this point of view, we expect that the 
same phenomena may be observed when the v component of the zero state is pointwise stable and not just 
convectively stable. 

We give the following example, 

u t = du xx + a (u — m 2 ) + fiv(l — u) 

Vt = v xx -^v. (4.1) 

We consider compactly supported, positive perturbations of the homogeneous zero state. With this initial 
data, v(t,x) will decay to zero for all x as t — > oo. With v converging to zero, one might expect that the 
u component will then spread with the speed 2y/da. However, this is not the case as one can observe in 
Figure [3j The decay rates of v(t,x) as x — > ±00 play an important role in the dynamics despite the fact 
that the solution is very small and converging to zero. This should be contrasted with the dynamics where 
v is set identically equal to zero. Here the u component converges, as expected, to the Fisher-KPP front 
traveling with speed 2yda. 



Figure 3: Anomalous spreading for (4.1) with d — 0.001, a 



= 4, = 2 and 7 = 0.1. The solid line is 
Note that the v component relaxes to zero 
0, 6, 12, 18, 24, 30, 36. On the left the 



the u component while the dashed line is the v component, 
while the u component spreads. The evolution is shown at times t 
evolution is shown for the full system. On the right, the same simulation is run but with identically zero 
initial conditions for the v component. The spreading speed in the absence of coupling is 2V0.004 s» .1265 
whereas the anomalous speed is approximately 1.9765. 



This is completely analogous to the unstable case described above. The linear spreading speed can be 
determined from the dispersion relation, as before, however the spatial eigenvalues for the v component are 
now 



± 



+ 47 + 4A. 



2 2 

The evolution of the linearized dynamics are again prescribed by a pointwise Green's function and computing 
pinched double roots leads to a precise characterization of the linear spreading speed. In particular, we find 
that poles of the pointwise Green's function arise whenever v~ = 1/+ and for the wavespeed 




This anomalous speed exists for d < 1/2 and 



a > 



d:-f 



1 - 2d 



We remark that for (4.1), if a pinched double root occurs then it is of the appropriate form to lead to 



faster speeds of propagation in the nonlinear system. If the coupling were in the opposite direction, that is if 
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the P term were in the v equation rather than the u equation then this pinched root would once again exist 
and give rise to an anomalous spreading speed in the linear equation, but would not lead to faster speeds of 
propagation in the nonlinear system. This is what occurs in ( 1.2 ) and explains why the fast linear spreading 
that was observed in that model does not persist to the nonlinear regime, see |13) . 



5 Numerical simulations 

In this final section, we compare numerically derived spreading speeds for the full nonlinear system to the 
linear spreading speeds predicted in Theorem [T] Recalling our discussion in section [3j we expect to observe 
anomalous nonlinear spreading only when d < 1 and within the parameter regime detailed in Theorem [T] 
We compute spreading speeds using a Crank-Nicolson type finite difference method. Initial data was chosen 
to be step functions in both the u component and the v component. Altering the initial location of these 
step functions relative to one another does not influence the selected wavespeed. Numerically observed 
wavespeeds for the u component are presented in Figure [4] and Figure [5j 




Figure 4: On the left, we plot numerically observed spreading speeds for the u component for d € [0, 3] and 
a 6 [0,3]. Note that the spreading speed remains two in the linearly anomalous regime for d > 1. On the 
right, we plot the relative error s b se rved/ max{2, 2\/da}. 



We remark on two qualitative observations. Throughout this article we have considered exclusively the 
case of f3 — 2. It is interesting to note that the particular value of j3 > has no influence on the wavespeed 
selection. Instead, we see that /3 influences which particular translate of the traveling front that is observed. 
Furthermore, this translate of the traveling front does not depend on the initial location of the u step 
function. The selected u front is determined solely by the v component and, after some initial transient, any 
compactly supported initial data for the u component will converge to this translate. 




Figure 5: Here we compare numerically derived spreading speeds within two separate sections of the (d, a) 
parameter plane. On the left, we consider a = 1.75 and vary d between and 1.4. The * are numerically 
derived spreading speeds while the curved dashed line is the anomalous speed, the solid line is the speed 
s u ,Un and the straght dashed line is s Vt u n . On the right, we do the same comparison, but now with d = 1.75 
and a varied. As predicted, these parameters lie in the region where there is no anomalous spreading. 
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6 Discussion 



In this article, we have investigated anomalous spreading speeds in a system of coupled Fisher-KPP equa- 
tions. We conclude with a brief commentary on more general aspects of this phenomena as well as a discussion 
of unresolved issues. 

The fact that the v component completely decouples in the nonlinear system is not a prerequisite 
for anomalous spreading. For example, the nonlinear system 

Ut = du xx + a (u — u 2 ) + flv(l — u) 
Vt = v xx + euv(l - u) + (y - v 2 ) , 

exhibits nonlinear anomalous spreading where the v component travels at speed 2 and the u component 
travels at the speed selected in the decoupled (e = 0) case. We note that qualitatively similar invasion 
phenomena has been observed in [2] in the context of an ecological model of the interaction of healthy and 
diseased red and grey squirrels. This model is much more complicated than ( |1.1[ ), but numerical simulations 
reveal behavior reminiscent of anomalous spreading. 

One might also be interested in the behavior of the system under some small perturbation that destroys 
the skew-product nature of the linearization. For example, consider the system 

u t = du xx + a (u — u 2 ) + f3v(l — u) 

Vt = v xx + eu(l - u) + (v - v 2 ) . (6.1) 



When e = 0, we recover (1.1 1. However, for e ^ 0, the linearization about the unstable state no longer 
has the requisite skew-product structure. In this case, the pinched double root of the dispersion relation 
that led to a pole of the pointwise Green's function perturbs to a branch point wherein G\ loses analyticity. 
These branch points are then dynamically relevant irregardless of whether they corresponded to nonlinearly 
relevant anomalous spreading speeds in the e = case. That is, for e small and positive the observed 



spreading speeds for (6.1) are small perturbations of the linear spreading speeds given in Theorem [T] and 



plotted in the left panel of Figure [T] If e < 0, then positivity is not preserved for (6.1 ) and solutions diverge 



However, if the quadratic nonlinear terms are replaced with cubic ones then for e < the linear spreading 
speed corresponds to values of A* for which Im(X*) 0. In this regime, the invasion front forms patterns 
in its wake with long wavelengths. See [10 for a precise example of this phenomena in the context of a 



phase- field model. Thus, system (1.1 1 sits on the boundary between systems with invasion fronts that leave 
patterns in their wake and those that don't. 

To conclude, we remark that a rigorous mathematical description of the anomalous invasion process 
remains an open question. As we pointed out above, the use of a traveling front approach is complicated by 
several factors. Most glaringly, is that the dynamics of the full system evolve on two separate time scales and 
therefore no choice of the traveling wave coordinate reduces the problem to an ordinary differential equation. 
More subtle, is the fact that the observed spreading speed depends on the "transient" evolution of compactly 
supported initial data towards a traveling front profile. 
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